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L p - VARIATIONS FOR MULTIFRACTAL FRACTIONAL RANDOM 

WALKS 

By Carenne Ludena 

Institute) Venezolano de Investigaciones Cientificas 

A multifractal random walk (MRW) is denned by a Brownian mo- 
tion subordinated by a class of continuous multifractal random mea- 
sures M[0,t], < t < 1. In this paper we obtain an extension of this 
process, referred to as multifractal fractional random walk (MFRW) , 
by considering the limit in distribution of a sequence of conditionally 
Gaussian processes. These conditional processes are defined as inte- 
grals with respect to fractional Brownian motion and convergence 
is seen to hold under certain conditions relating the self-similarity 
(Hurst) exponent of the fBm to the parameters defining the multi- 
fractal random measure M. As a result, a larger class of models is ob- 
tained, whose fine scale (scaling) structure is then analyzed in terms 
of the empirical structure functions. Implications for the analysis and 
inference of multifractal exponents from data, namely, confidence in- 
tervals, are also provided. 

1. Introduction. Multifractal models, or more generally, multifractal ran- 
dom walks, with respect to ordinary Brownian motion or fractional Brown- 
ian motion have become very popular models for scaling phenomena. This 
kind of phenomena is associated to very nonlinear models as encountered, 
for example, in hydrodynamic turbulence, biological rhythms or the network 
traffic [4, 5, 20, 21]. 

The scaling property is usually defined in terms of a scaling exponent, 
Cip), associated to the process X(t) describing the phenomena. This expo- 
nent is defined as 

_ logE(\X(t + r)-X(t)\P) 
^ P) - log(r) 
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at different scales < r < T. Based on the multifractal formalism [2, 14, 19] 
which establishes the relationship between £(p) and the characterization 
of the fluctuations of the local regularity of X(t), estimating this scaling 
function has attracted a great deal of attention both from a theoretical and 
a practical point of view (see, e.g., [1, 16, 23]). 

Typically, given a sample path of X(t), t G [0, 1], say, estimators of Q{p) 
are based on the logarithm of the empirical structure function of the process 
over a certain regular partition to, . . . , t m of [0, 1] (ij — tj_i = t), that is, of 



fc=i 

This procedure implicitly assumes that S m {jp) converges a.s. to E(|X(t + r) — 
X(t)\ p ) as m->oo (i.e., ergodicity), and thus, that log(5 m (p))/log(r) is a 
consistent estimator of C(p). We will refer to this method as a moment or 
L p -variation estimator. 

L p -variation estimators have two major drawbacks in practice. On the 
one hand, they tend to be very data demanding [23]. On the other hand, 
out of a certain range [p_,p+] a certain linearization effect ([16, 22, 23] 
and the references therein) tends to appear. That is, log5 m (p)/log(r) — ► 
kp for a certain constant k instead of the expected C(p)- This has led to 
consider estimators based not on the increments, but rather on the suprema 
of the increments or likewise techniques which provide more robust methods 
[2, 16, 23], although this phenomena is not thoroughly understood. 

The aim of this paper is twofold. In the first place, based on the con- 
struction of multifractal random measures (MRM) developed by Bacry and 
Muzy in [7] , we consider an extension of their multifractal random walks for 
ordinary Brownian motion to multifractal fractional random walks (MFRW) 
based on fractional Brownian motion with Hurst parameter H. Two possible 
courses of action are proposed. In each case the existence of the process is 
assured under certain additional conditions involving the Hurst parameter 
and the multifractal random measure. 

In the second place, we study the conditional distribution, given the mul- 
tifractal random measure, of the structure function for even and positive p. 
As a corollary, we discuss the properties of a L p -variation estimator of the 
scaling function, for even p. The stated conditional asymptotic distribution 
allows us to construct confidence intervals for for large values of m, the 
number of observations. Our main result shows that approximation errors 
are fundamentally due to a bias term which tends to zero at a logarithmic 
rate. 

On the other hand, a very interesting phenomena, related to a certain 
linearization effect reported in many applications, appears for intermediate 
values of m. This phenomena is related to the range of scales over which the 




//-VARIATIONS FOR MFRW 



3 



scale invariance property holds true. In the MRM model developed in [7] this 
range depends on a certain parameter T. Based on this characterization, we 
show that the behavior of the estimator will be very different according to 
the relative size of T and m. If T <C m (T = 1, say), the asymptotic behavior 
of the estimator will depend on the conditional asymptotic distribution of 
the structure function. However, if T ~ m, for big enough m, the scaling 
property will introduce a linear bias term in the estimation procedure. 

The article is organized as follows. In Section 2 we introduce the class of 
MRM following the construction of Bacry and Muzy and define the MFRW 
based on fractional Brownian motion with Hurst parameter H > 1/2. In 
Section 3 we study the asymptotic properties of the sums of certain functions 
of the increments of the MFRW for 1/2 < H < 3/4. In Section 4 we apply the 
asymptotic results of Section 3 to study the properties of a moment based 
estimator of the scaling function. In Section 4.2 we discuss the problem of 
estimating £(p) when T ~ to, that is, when the number of observations is 
comparable to the scale invariance range. 

2. Multifractal random measures. The existence and analysis of discrete 
multifractal random measures, based on the Mandelbrot multiplicative cas- 
cades [17], was first established by Kahane and Peyriere [15]. Continuous ex- 
tensions were introduced by Barral and Mandelbrot [8] . This class of models 
was further extended by Bacry and Muzy [7] who define a class of stationary 
random measures (MRM) based on the following elements: 

1. Let v be a Levy measure defined over IR satisfying 



3. Define the process P(A),A C S + , such that, for any sequence of dis- 
joint sets A n C S + , (P(A n ))^ =1 are independent r.v., POJ^An) = 




for all u > 0. Associated to v and m£K, define for q £ M 



, , N f e q — 1 — iq sin x . , . 

(p{q)=imq+ / ^ u{dx). 

J x l 



2. Set S + = {(t,l),t <EM + } and the associated measure 

fi(dt,dl) =r 2 dtdl. 



E~=i P{An) a.s. and 



e <t>{qMA)_ 



4. Let 



/(*) 




T. 



k<T 
k>T 
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where T is arbitrarily large. For t E R and I > 0, define the collection of 
sets 

(1) Mt) = {(s, k), k>l, -f(k)/2 <s-t< f(k)/2}. 

Next for iGR, consider the process wi(t) = P(Ai(t)). 
5. For / > and / a Lebesgue measurable set, define 

M^I) = ^e w ' {t) dt, 

and set Mi(dt) to be the associated measure. It can be seen (see exact 
conditions and references below) that there exists a random measure 
M(dt) such that a.e. Mi(dt) converges weakly toward M(dt). The MRM 
measure is then defined as the limit 

M{dt) = lim Mi(dt). 

The properties, and existence, of M(dt) are discussed in [7], as well as several 
examples of popular choices of the measure v. 

Associated to the set process P and <fi, we may also consider the real func- 
tion ip(q) = <j)(—iq) (see [7] for further details) defined in terms of the moment 
generating function, E(e qP ^) = e ^^ A \ for any AcS + . Throughout the 
article we shall assume = 0. Let ((q) = q — ?p(q), and introduce the 
following assumptions: 

[Al] There exists e > such that ((1 + s) > 1 
[A(?)] C(?)>1- 

The following lemma provides a summary of the properties of MRM ob- 
tained in [7] and which we list here for ease of reference. 

Lemma 2.1 (Theorems 2, 3, 4 and Lemma 5 in [7]). 

1. [Nondegeneracy of M(dt).J If there exists e > such that £(1 + e) > 1 
([Al]), then E(M([0,t]))=t. 

2. (Positive order moments of M .) If £(g) > 1 ([A(q)]) for q > 0, then 
E(M([0,t]) q ) < oo. Moreover, if [Al] holds and E(M([0,t]) q ) < oo, then 

C(?)>1- 

3. [Moments of sup ue [ t ] e Wl( ~ u \] If ip(q) ^ oo [in particular, if ((q) > 1], 
then 

e( sup e qw '^A <oo. 
\ue[o,t] / 
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4. (Scaling property.) For all t <T, 




whenever £(g) > 1 . Moreover, 



{wix(Xt)} t 



law 



n x + { Wl {t)} t 



and 



{M[0,Xt}} t ^ \e^{M[0,t}} t 



for all A G (0, 1), where Q\ is an infinitely divisible random variable which 
is independent of {wi(t)}ij with characteristic function 

E(e iqnx ) = X'^. 

2.1. Constructing multifractal fractional random walks. Assume through- 
out this subsection, without loss of generality and in order to simplify the 
construction of the MFRW, that T = 1. For M a nondegenerate measure 
(i.e., such that E(M([0,i])) = t), Bacry and Muzy [7] define the subordinate 
random walk as 



for B(t) a Brownian motion which is independent of M. This process can 
be interpreted as the centered, conditionally Gaussian process whose condi- 
tional covariance function is defined as the limit 



The authors also show this process is equivalent in law to the limit process 



This procedure has been repeated based on heuristic arguments for frac- 
tional Brownian motion (see, e.g., [6, 16]). The work by Bacry and Muzy 
[7] provides the basis for a rigorous definition based on the properties of 
the stochastic integral with respect to fBm whenever the integrand is de- 
terministic (see, e.g., [10]) and the Hurst coefficient, H, is greater than 1/2. 
However, though it might appear to be the natural thing to do, the MFRW 
cannot be defined directly as the limit in distribution of the (well-defined) 
stochastic integrals Jq e Wl ^ s ^ 2 dB^ . Indeed, although for each I > 0, the con- 
ditional covariance Ri(s,t) associated to the stochastic integral, with respect 
to fBm for H > 1/2, is well defined and given by 



Y t = B(M([0,t])) 



T(s,t) := lim / 
i^o+ Jo 





R l {s,t)=C H l [ 



s 



U — V 



2H-2 e Mu)/2 e Mv)/2 dudV) 



Jo Jo 
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a direct calculation of the expectation using the moment generating function 
of process {wi}t (Lemma 1 in [7]), yields that liminf; KRi(s, t) — > 0. More 
precisely (assume s < t), we have 



E J R z ( S ,t) = e 2 ^ 1 /2) r 2^(i/2) 



t rs 



JO 



u — V 



\2H 



-2 e -2^l/2) Pl {u- V ) dudv 



2H 



V 



-2 e -2f{l/2) Pl (v) dvdu 



JO 

t rit 



+ 



; 2H-2 e -2^(l/2) ft («) dvdu )/ 1^(1/2) 



r 2^(1/2) 2 



I ru 



JO 



v~~ -e-WO-mPiMdvdu 



.2H-2, 



+ 



I ru 

+ / 

Jl 



t ru 



where 



+ ^(1/2)^-2^(1/2) 

I + II + III + IV, 



( ln(l/Z) + 1 - v/l, 
Pi{v) = l In(l/u), 



v 2H-2 e -2i,{l/2) Pl {v) dvdu 



v 2H-2 e -2i,(l/2) Pl (v) dvdu 



10, 



if v < I, 
if l<v<l, 
if v > 1. 



Setting = 2-0(1/2) 



/ = 2Z 2 " (Y i 2 ^ 2 (i - l)e at dt^j - 



and 



// = 2(s - Z)^ -1 / X t 2H_2 (t - l)e a * dt 
On the other hand, for Z < « < 1, we have p;(w) = log(l/u), so that 



///<2e a (r 2 ^ 1 / 2 ) / t 2H-2 + 2^l/2)+l dt _i2H-l 







since H > 1/2 and ■0(1/2) < by convexity. Finally, /V can be bounded as 
above to show it also tends to zero. 

Set R(s,t) := lim^ + Ri(s,t). An application of Fatou's lemma yields 
ER(s, t) = 0, and whence that R(s, t) = 0. 

In what follows we propose two alternative procedures. In each case an ad- 
ditional condition relating H and the scaling function ip(-) must be imposed 
in order to assure the existence of the limiting process. 
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2.1.1. A first approach. Start by defining the measure MRM M 1 / 2 , when- 
ever the scaling function ip(l/2) > —1/2, as follows. For any Lebesgue mea- 
surable set /, M X ' 2 (I) = lim^ - Z^(V2) e -iKt)/2 e w,(*)/2 dL Existence of this 
measure follows as for the proof of Theorem 1 in [7] . 

With this construction, we have the following definition of a MFRW, 
where here and below Em denotes the conditional expectation given the 
MRM M. 

Definition 2.1. Let H > 1/2. The MFRW-1/2 is the centered process 
Xt which, conditionally to M(dt), has a Gaussian distribution with condi- 
tional covariance 

S{s,t) = C H C ( S \u-v\ 2H - 2 M 1 l 2 {du)M l l 2 {dv), 
Jo Jo 

where Ch is a constant which only depends on H. 

This conditional covariance can be obtained as the limit S(s,t) = 
linn_ >0 + S[(s,t), where 

(2) S^-Ch f [ S \u-v\ 2H - 2 l- 2 ^e w ^l 2 e w ^l 2 dudv. 
Jo Jo 

We have the following lemma, which assures this limiting covariance is 
well defined if (Hl/2): H + ip (1/2) > 1/2 holds. Remark that because tp(-) 
is convex and if)(0) = ^(1) = 0, then ip(l/2) < 0. 

Lemma 2.2. Assume (Hl/2) holds. Then S(s,t) exists a.e. and is not 
identically equal to zero. 

Proof. Using the moment generating function of process {wi}t (Lemma 
1 in [7]), for u<v, 

E( e l/2(>i (u)+wi (y))\ = e 2^(l) Pi (0)+(V(l)-2^(l/2))p,(^-«) 5 

where, by construction, ^(1) = 0. As pi(0) = log(l//) + 1, we have 

E(S l (s,t)) = C H f /%-^-2 e -2^(i/2)MI«--l). 
jo Jo 

On the other hand, pi(u) = log(l/n) + log(ti//) + 1 — u/l for < u < I, and 
pi(u) = log(l/n) for I < u < 1, so that, for u > 0, 

e -2f(l/2) Pl (u) < e -2V(l/2)log(l/ M ) 

for all / > 0. An application of the Dominated Convergence Theorem yields 
that the limiting random covariance function exists a.e. 
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For the second part of the lemma, by direct calculation of the expectation 
as above, sup^ E(M ; 1/2 ([0, t])) < oo and E(M 1 / 2 ([0, t])) < oo since ^(1/2) > 
— 1/2 as follows from (Hl/2). As in the proof of Theorem 2 (Lemma 3) in [7], 
this yields the limiting measure M 1 / 2 is nondegenerate (E(M 1 / 2 ([0,t])) = t). 
For any < a < s A t and < e < s At — a, 

S(s,t)> P f lu-^^M^^M 1 / 2 ^) 

> £ 2H-2 M l/2 ([0)a])M l/2 ([a + £)S])) 

where the last expression is not zero a.e. □ 

Hence, the conditional covariance is well defined, as well as the process 

x t . 

On the other hand, let B H (t) be a fractional Brownian motion, with 
H > 1/2 and E(B H (1)) = 1, which is independent of the MRM M. It is also 
possible to define the MFRW-1/2 process directly, as the limit of 

X lt := / t r^ 1 /2) e ««(u)/2 dB H ) 

Jo 

as I — > + . This integral is understood, conditional to M, in terms of the 
Gaussian construction for deterministic integrands. That is, for each I > 0, 
X\j\m is a conditionally centered Gaussian process with covariance Si(s,t). 

Let X' t be the limiting process whenever it exists. The following result 
shows that both procedures are equivalent. 

Lemma 2.3. Assume (Hl/2) holds, then X t = X[. 

Proof. That the finite dimensional distributions of X\ t converge to 
those of Xt follows directly from the definition so that in order to show the 
stated result it is enough to see that the sequence is tight in the Skorohod 
metric. For this, we start by observing that 

E(|X Ml -X M2 | 2 ) 

= E(C H \U - v fH-2 r 2f(l/2) e Muy2 e m(v)/2 dudy 

as follows by calculating the conditional expectation. 
Hence, 

(3) E(\X lM -X ljt2 \ 2 ) 

(4) = C H f 2 \U - v fH-2 e -2W/2) Pl (\u-v\) dudv 
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<C H [ 2 \u-v\ 2H - 2+2 ^dudv 

(5) =C(H)\t 1 -t 2 \ 2 ( H - 1+ ^ 2 » +2 , 

for a certain constant C(H) and t2,t± £ [0, 1], which proves tightness under 
the stated condition. □ 

2.1.2. A second approach. 

Definition 2.2. Let H > 1/2. The MFRW-1 is the centered process X t 
which, conditionally to M(dt), has a Gaussian distribution with conditional 
covariance 

R{s,t) = C H f I" \u-v\ 2H - 2 M(du)M{dv), 
Jo Jo 

where Cjj is a constant which only depends on H. 

This conditional covariance can be obtained as the limit R(s,t) = 
limj_ >0 + Ri(s,t), where 

(6) Ri{s,t):=C H f I" \ u - v \ 2H - 2 e w ^ u) e w ^ v) dudv. 

Jo Jo 

The existence of this limiting covariance is assured if the following condi- 
tion (HI): H - V>(2)/2 > 1/2 holds. Convexity of </>(•) yields ip(2) > 0. 

Lemma 2.4. Assume HI holds. Then R(s,t) exists a.e. and is not iden- 
tically equal to zero. 

Proof. The proof follows exactly as that of Lemma 2.2. Nondegeneracy 
of M(dt) follows from HI, as H - ip(2)/2 > 1/2 yields £(2) > 1. □ 

Hence, the conditional covariance is well defined, as well as the process 
X t . 

On the other hand, let B H (t) be a fractional Brownian motion, with 
H > 1/2 and E(B H (1)) = 1, which is independent of the MRM M. It is also 
possible to define the MFRW-1 process directly, as the limit of 

X ltt := fe w ^dB* : 
Jo 

as I — > + . This integral is understood, conditional to M, in terms of the 
Gaussian construction for deterministic integrands. That is, for each I > 0, 
Xij\m is a conditionally centered Gaussian process with covariance 

R l (s,t) = C H f I" \u - v\ 2H - 2 e w ' {u) e Wl( - v) dudv. 
Jo Jo 



10 



C. LUDENA 



Let X' t be the limiting process whenever it exists. The following result 
shows that both procedures are equivalent. 

Lemma 2.5. Assume HI holds, then X t = X[. 

Proof. The proof follows exactly as that of Lemma 2.3. □ 

2.1.3. Existence of moments. An important issue is stating sufficient 
conditions for the existence of the moments of the MFRW-1/2 or the MFRW- 
1 process Xt- In each case these will be related to the scaling function ?/>(•). 

Introduce the following definitions: 

Definition 2.3. We will say process satisfies assumption Hl/2-p if the 
MFRW-1/2 process X t satisfies 

E(|X 1 -X | p )<oo. 

Definition 2.4. We will say process satisfies assumption Hl-p if the 
MFRW-1 process Xt satisfies 

E|Jfi - X \ p < oo. 

Although both constructions are roughly equivalent, the scaling factor will 
not be the same in every case. Indeed, let Xt be a MFRW-1/2. If 1 < j < n, 
the scaling property (4, in Lemma 2.1) yields 

E(|X (j _ 1)/n - X j/n \ 2p ) = ^W-W/2)-2pff E( | Xi _ ^py 

If, however, Xt is a MFRW-1 process, then we obtain the more familiar 

E(|X _ 1)/n - X l/n \ 2p ) = n^~ 2pH n\Xi - X \ 2p )- 

In what follows we will restrict our attention to the second approach. 

Set 5 m := ip{m) + ip(m — 2) — 2ip{m — 1) for m > 2. Convexity of function 
ijj assures that 5 m > for each m > 2. Remark that in the case m = 2, 
^2 = '0(2)- The next lemma gives sufficient conditions for Hl-2p. Its proof 
is based on power counting techniques, namely, Section 3 (Theorem 3.1) in 
[12]. Power counting techniques provide sufficient conditions for the existence 
of integrals whose integrands are products of powers of affine functionals. 
Based on these results we have the following lemma. 

Lemma 2.6. Assume 2H — 1 — S m > for each 2 < m < 2p, then as- 
sumption Hl-2p holds true for p6N. 
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Proof. Calculating the conditional expectation, we have, for p £ N, 
E(|AV - X lfi \ 2 P) = cp&(ff[ Wi ~ v^^M+mM dui . . . dv \ . 

Consider any arbitrary permutation of {ui , v\ , . . . , u p , v p } given by {z\ , . . . , 
Z2p}. As follows from the properties of the moment generating function of 
wi, over the set z\ < ■ • ■ < z 2p the last integral is bounded by 

/ n n i* - Zj \ h ^ 2H -v-^ tfai • ■ • d Z2p , 

J i=2j<i 

where h(i,j) 6 {0,1} and h(i,j) = 1 exactly p times according to whether 
Zi = Uk and Zj = in the original order [of all the p(2p — 1) couples there 
are exactly p that satisfy this condition]. Following the notation of Section 
3 in [12], we are interested in the p(2p — 1) linear functionals Lij = (zi — Zj) 
with j <i. 

By the power counting theorems, existence of the integral will follow if 
for any strongly independent subset of size k of these linear functionals the 
condition 

k 

k + J2 h dhji)(2H - 2) + 5 k - jl+1 > 
i=i 

holds true. In this setting any strongly independent subset of linear function- 
als can be identified as a tree (or a collection of disjoint trees) constructed 
over the set of vertex {zi, . . . , z 2p } identifying each pair {zi,Zj} with an 
edge. It follows that the maximum size of any strongly independent subset 
is 2p — 1 . 

If k < p and if we choose the k terms with h(i[,ji) = 1, then 

k k 

k + £ h(i h ji)(2H - 2) - 6 k - jl+1 = J2l(2H - 1) - > 

i=i i=i 

under the assumption that (2H — 1) — 5 m > for each 2 < m < 2p. 

If k > p or if we choose a subset with h(ii,j{) = 0, the result follows as 
1 - S m > 2H - 1 - 6 m . □ 

Remark 2.1. When H = 1/2, Hl-2p is just A2p and is assured whenever 
2p - 1 - tp(2p) > 0. If p = 2, Hl-2 is assured if 2H - 1 - ip(2) > 0. It is 
hence fair to conjecture that a sufficient condition for Hl-2p is given by 
2pH — 1 — ip{2p) > 0, which in turn assures together with 2H — 1 — ip{2) > 0, 
that kH — 1 — tp(k) > for any 2 < k < 2p, as follows by the properties of 
function if). 
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It is important to remark that in any case Xt is not the subordinated pro- 
cess U t = B H (M[0,t]). The latter indeed is a centered, conditionally Gaus- 
sian process, with conditional covariance 

Cov(C/ t ,[/ s ) = i[(M[0,t]) 2 ^ + (M[0, S ]) 2H -|M[0,t]-M[0, S ]| 2H ]. 

This process is referred to as fractional Brownian motion in multifractal 
time [13, 16]. 

3. Quadratic and higher-order moments of the increments. Formally, 
we may be tempted to assume dXt = lim;^ + e Wl ^dB^ , however, u>i is not 
necessarily regular and this interpretation may lead to confusion if trying 
to interpret the above expression by considering pathwise integrals. We will 
thus consider Xt in terms of the conditionally Gaussian specification. In 
particular, throughout the rest of the paper we will assume HI holds true. 

Throughout the next sections we will consider dyadic increments, so al- 
ways m n = 2 n . We also assume, in order to simplify notation, that T = 1. 
Assume we observe X tl , . . . , X tmn with tj = j/m n . As mentioned in Section 
1, usual procedure for estimating is to look at the average of the p 
moments of the increments (or, more generally, the increments of order r, 
which are roughly equivalent to considering wavelet coefficients for a wavelet 
with r vanishing moments) and then use the scaling properties (see [1]). Our 
main interest in this section will be developing an insight as to the estimation 
problem at hand, by looking at the asymptotic behavior of these averages for 
even p and r = 1. Because of the long range dependence property of fBm for 
H > 1/2, it is necessary to further restrict the value of the Hurst parameter 
by asking that H < 3/4 in order to obtain a limiting conditional Gaussian 
behavior. This restriction can be avoided by considering the r increments 
of X(t), with r > 1. However, this involves further technicalities which were 
not the main interest of this work. 

We start by defining the increments of Xt. For j = 1, . . . , m n , define AXj = 
X j/m n ~ x (j-i)/m n with Xq = 0. 

Set dj^n = EA/[(AXj) 2 ] 1 ^ 2 , and define 



Z n (p) = — £[|AX,|* '- c^J, 



'i=i 



where c p = E(|Z| P ), Z ~ N(0, 1). 

Finally, define the conditional correlation 

(v -\ f A \" a y i ^ijAXjAX, 



k 



0'j,nO'k,n 

We are interested in studying the asymptotic behavior of Z n (p), based on 
the scaling property and the conditionally Gaussian structure. 
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By construction, 



u-v\ 2H - 2 M(du)M(dv) 



and 



-v\™-*M(du)M{dv) 



1/2 



2H-2 , 



Q , j,nQ j k,n 

As mentioned above, the asymptotic behavior of Z n relies on the condi- 
tionally Gaussian behavior of process Xf. The next two technical lemmas 
provide the necessary tools. 

Lemma 3.1. There exists a (deterministic) constant C, such that 
(7) p(j,k)<C\j-k\ 2H - 2 , 
for j / k. 

Proof. Assume to begin with that \j — k\ > 2. Remark in p(j,k), the 
numerator is bounded by 



k-j-l\ 2H ~ 2 

n / 



m. 



tj-i 



M(du) / M(du) 



On the other hand, we may bound ai n from below 



(8) 



1 

^ > ~m m ~- 
1 1 in 



ti-i 



M(du), 



for 1 < I < m n , which yields the stated bound. For \j — k\ = 1, we have 
p(j, k) < 1 since p is a correlation. □ 



Define for r E N, r > 2 

r„(r,p) 



in, 



J2 P( k >3Y a j,n a k 



k,n' 



l<fc,j<m„ 



Lemma 3.2. Assume 1/2 < H < 3/4 andp E N. T/ien ; m„ p VWP T n (r,p) 
converges a.s. to a positive integrable r.v. T(r,p) for all even r <p, whenever 
p is even, Hl-2p holds and condition A(p): 2pip'(2p) < ip(2p) + 1 is satisfied. 

PROOF. By stationarity, E(a? n a£ n ) <E(al P n ). And as follows by the 

scaling property, this last expectation is equal to mt, 2p ^E(|Xi — Ao| 2p ), 
which justifies the normalizing constant. 
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On the other hand, Lemma 3.1 guarantees that p(k,j) r /C r \j — k\ T ( 2H ~ 2 ^ < 

1. 

Set a_j jn = a,j >n and define, for \u\ < m n , 

w ^ ) m ^' P (k,k + uY p p 

Wu ' n ~ m n -\u\ ^ max(\uliy(^-2) a k,n a k+u,n- 

Hence, 

m^^)r n (r,p)= £ (l-MW(^-2) ]w ^ 

|«|<mn-l V 

We will start by showing that W U)Tl converges a.s. to a certain integrable r.v. 
W Uj r for each fixed u, \u\ < m n . 

I. As part of the proof, we will require a truncation argument which holds 
for all \u\ < m n simultaneously. Let r > 1 be such that a' := ip(2pr)/r — 
tp(2p) < y/2 — 1/r — 1. This can be seen to hold for r close enough to 1, 
by checking that F(l) = and F'(l) < for F(r) = i(j(2pr) - ri/j(2p) - 
ry2 — 1/r + r under A(p). Since \J2 — 1/r — 1 < 1 — 1/r, this choice of a' 
yields that a := i^(2pr)/r - ip(2p) + 1/r < 1. Set 7 = (1/4 + 3a/4), so that 
a < 7 < 1. Set Aj = {ml pH -^ 2p) a^ < ml} and A n = Aj. 

Then, by the scaling property, 



m " E(m r n {2pH -^ 2p)) a^) 



<r" /m 1 + 1 /'( 2 P T ') — f^(2p)— r7 _ nnr(a— 1)/4 
— ' u n ~~ z J 

which is the term of a convergent series. Thus, by the Borel-Cantelli lemma, 
P(A n i.o.) = 0. 

II. To show a.s. convergence of W UtTl , we start with u = 0. Remark in this 
case max (Q^r(2H-2) = 1- The proof follows from a series of steps. The ideas 
behind this part of the proof are closely related to the proof of Theorem 3.1 
in [18]. 

Step 1. Set l n = 2-( 1 " a ) n , where a = (1 - o)/2, and let .F n = a{w s (u),s > 
l n , t £ [0, 1]}. We shall see E^ n (W ,n) - W , n -> a.s. 

Set Aj >n = [(j — l)/m n , j/m n ], let 2; = (iti,wi, . . • , u p , v p ) and define, for 
any < I, 



2pH-i/j(2p) p 

<m*)=^ — En 



m ™ V TT l Ui - Ui | 2 ^- 2 e u ' i ^)- ,i;i «^)e u ' i(,;i) "' Ui " (,;i:) 



J 2=1 



/ p \ p 

E [ JJ e «>«n(«i) e ««n(«i) m^A 2 (Ui,Vi) 

\i=i J i=i 3 ' n 
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2pH-i,{2p) p 

e n ,i{z) = y,\\ \v>i-Vi\ 1 A 2 {ui,Vi) 

n 3 *=i 



Mil 



\i=l 



It is important to remark that if I < l n , for each n, wi(s) — wi n (s) is 
independent of wi n (t), for each s,t £ [0,1], because, by construction, the 
sets are disjoint. Hence, 



e n 



n 

i=i 



(ui) e Wi(Vi) 



: ]W TJ e »l n ("i) e »l n ("i) W j TT e «'l(lii)-W In (tti) e 1fl l (« i )-«J In (« i ) 

\i=l / \i=l y 



and e nj ;(z) = E(d ni z(z)) so that M(J A d n: i(z)) = J A e n j(z) for any Borel set A. 
Moreover, let eo(p) :=lim^_ >0 + / e n,l( z )dz, so that 

e (p)=E^jf[\u i -v t \ 2H - 2 dM(z) ) j, 

as follows by the scaling property and Lemma 2.5 whenever HI holds. 

On the other hand, because the sequence m n is dyadic, e Wl ^ u ' Wl ^ u > and 
e v)i(v)-v>i n (v) are independent if / < l n and u E Aj jn and u G Afc jn with |j — 
fc| >2 na . 

Set 

TIP p wi n (u,)+w ln {vi) 

h n (z):- lU=1 



E(nf=l e^n^+Wln^i)) ' 



With the above notation, note that Wo, n = hrm_ >0 + / h n {z)d n j(z) dz and 
Ejr n (W 0in ) = lim^ + / h n (z)e ny i{z) dz. Now consider 



A n := lim / h n (z)(d n j(z) - e n j(z)) dz 



and 



A n :=l A ^lim / K(z)(d n< i(z) - e n j(z)) dz. 
First, since P(A„ / A n ) = P(A^), it is enough by I to study A n . 
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Second, by calculating conditional expectations and recalling the terms 
in the sum are conditionally 2 na -dependent, we have 

E(A 2 n ) < E ™r^ (2p) E(af ) = 0(mr 1+a ), 

mi. ^ J 

n J= i 

which is the term of a convergent series in n by the choice of 7 and a. 

Step 2. We shall now prove that lim;_ >0 + / h n (z)e n j(z) dz converges a.s. 

First remark that, for each fixed I, wi(u) — wi(0) = P(Ai(u) — -Az(O)) — 
P(Ai(0) - Ai(u)), where A - B = An B c for any two sets A and B. By 
direct calculation, recalling T = 1, 



fi(dl,ds)= / fi(dl,ds) = u/l + u 

Mu)-M0) JAM-Mu) 

and 

Hence, by stationarity, if E(|^ n (0)| 9 ) < 00, there exists C > such that, 
if ue A J>5 then E(\w ln (u) - wi n (j / m n )\ q ) < C(u- j/m n )/l n < Cm~ a . 
On the other hand, calculating the expectation and the limit as I — > + , 

lim / h n (z)e n i(z) dz 

2 P H-i>(2p) 
\ I llijn 

^ m n 

m n p 

xY,J hn(z)Y[ 1 A? n { Ui ,Vi)\ui - Vi\ 2H ~ 2 Y[\zi - z k \~ a ^ dz, 

j=l i=l i,k 

where Yl' 1S the sum over all the possible (2p)\ orders of the 2p integrands. 
Let s stand for a fixed permutation, such that s(zi), s(z k ), i, k = 1, . . . ,p, are 
the given positions in the order. Then a(i, k) = ip(\s(zi) — s ( z k)\) +'4>(\s(zi) — 
s(z k )\+2)-2^{\s(z l )-s{z k )\ + l) and £ i>fc <*(», k) = ^(2p) (see [7]). 
Consider the expression 



2pH-^(2p) m n 

jj n ,-Y]'— V e 2 i™i„(iK) 

nf=i 1a? (ui,Vi)\ui - Vi\ 2H - 2 lii lk W - z k \~ a{i ' k) 



where, for each j, 

v 



Ui,Vi)\Ui - Vi\ 2H 2 



L p - VARIATIONS FOR MFRW 17 

*Y[\zi-z k \- a W dz = e (p). 
i . k 

And since there exists c > such that (Lemma 1 in [7]) 
E(U P =1 e w 'nM+ w 'n( v i)) > cm { n~ a)i,(2p) , the expression 



e'o(p) :=E 



t m 2pH-a^(2p) 



UU 1a2 {ui,Vi)\ui-Vi\ 2H ~ 2 n ilfc \zi-z k \ -<"&*) 

^— <C oc 

E(nf=i e^ ! " ( Mi ) +U, in (^) ) 

and does not depend on n, so that we can write 

Jo - 

and by Proposition 2.2 in [18], under A(p), U n converges a.s. to a certain 
random measure eo(p)M^[0, 1]. 

In order to end the proof of this step, set D n = linn__ >0 + / h n {z)e n j{z) dz — 
U n , which can be written as 

T-^ / ml pH ~' p( - 2p) ^ f e 2 P»"inU/mn) _ jj p =1 e ™i„(«i) 



nl A 2 (Ui,Vi)\Ui - Vi\ 2H ~ 2 Y[\zi - z k \ 

Set 



x 



=1 i,k 



(k-l)wi n (j/m n ) r-r 2 P p wi n 

Using the identity z\ . . . z p — w\ . . .w p = (z\ — W\)(z2 ■ ■ ■ z m ) + W\ [zi . . . z p 
W2 ■ ■ ■ w p ) and bounding 

<e^W|^„(«)-^ n (i/m n )|, 

we have 

2p-l 2pH~i,{2p) 

Ia^E'E — 



E y hiJ^l-^tiK) 



X 

3=1 
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p 



xe v>l n U/™ n ) b fc (s)TTi a ( Uu 
'=1 J '" 

xY[\zi-z k \- a ^dz. 



Vi)\Ui -Vi\ 2H 2 



i,k 

Let r > 1 be as in the proof of step 1. Again, by Lemma 1, [7], 

E{e rWl " ij/mn) b r k:n {z)) = 0( m a-«)W2pr)-r^(2 P )))_ 
Let r' = r/(r — 1). Holder's inequality yields 

ftH-Mp) f E{\ Wln {z k )- Wl Sj/m n )\e^l^\, n {z)) 
v 



■m: 



xTTl A 2 (u^Vi^iii-Vil 21 * 2 T[\zi-z k \ a{ i< k) dz 

i=l i,k 

<m 2 r P H -^ J [E(e rWl ^ j/mn) b r k>n (z))] 1/r 

^n\ Wln {zk)-w ln (j/m n )f)] l l r ' 
p 

x TT^A 2 {Ui,Vi)\Ui-Vi\ 2H - 2 

n\ Zi -z k \- a ^dz 



X 
i.k 

< Cm£~ a * > ^( 2pr ^~ r ^ 2p ^/ r m~ a l r ' 

— Qf m (l-<*)('ip(2pr)-r'ip(2p))/r-a(r-l)/r\ 

The exponent in the last expression can be rewritten as rj := (1 — a) a' — 
a{r — l)/r. In order to end the proof, it is enough to show that rj < and, 
thus, that the latter is the term of a convergent series. 

As follows from the choice of a, 

_ , (1-1/r-aQ 2 
V ~ a 2(a' + l/r) • 

Consider the second degree polynomial p(x) = 2x(x + 1/r) — (1 — 1/r — 
x) 2 = {x + y/2- l/r + l )(x + 1 - y/2 - 1/r). Hence, 77 = 2(a' + l/r)p(a') < 
since < a' < ^/2 — 1 Jr - 1 . 

III. Next, for u^O. For this part we require p to be even. We are interested 

2pH-4>(2p) 

in the convergence of mn J2kP r (^j^ + u ) a k a k+u- Since r and p are 

even, each of the terms in the sum can be written as the product of 2p 
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integrals 



p—r p—r 



al a 



k+u 



n 



Ui - V; 



\2H-2 



M(dui)M(dvi 



n 



\Ui - fj, 

n 

(p-r)/2 

>< n 



2H-2 



M(dui)M(dvi 



\Ui - Vi 



2H-2 



M(dui)M(dvi) 



i=i 

(p-r) 

n 

i=(p-r)/2+l UL ^k+u,n 



\Ui - Vi 



2H-2 



M(dui)M(dvi 



This product has exactly the same structure as that of the terms in Wq n , 
although, of course, the expectation will not be the same. 

On the other hand, applying the identity, 2ab < a 2 + b 2 , to the last product 
of integrals, 



2pH-^{2p) 



Y,P r (k,k + u)a p k a p k+u 



can be seen to be bounded over the set A n defined in I and the proof of step 
1 can be carried out analogously. For step 2, remark that 

E(m% H -^pr(k,k + u)ala p k+ J 

does not depend neither on k nor on n. The first assertion follows by station- 
arity. For the second, it is enough to consider a change of variables z[ = m n Zi. 
By Lemma 1 in [7], we have, for each fixed I > 0, 



2p 



Ki=l 



over the set u± < ■ ■ ■ < U2 P , where pi(u) = log(l/u) + log(u/Z) + 1 — u/l for 
< u < I and pi(u) = log(l/ii) for I < u < 1. Also, as J2i,j a (hj) = 
for each fixed / > 0, calculating expectations, we find a normalization factor 
m n 2pH+ ^'^ 2p ^ + ' 2p , so that by the change of variables, we obtain the factor 



m, 



2pH+i/i(2p) ^ jj ence calculating the limit expectation as / — > + , 
e u (r,p) := E(mf-*W/(fc, k + u)a p a p +u ) 



Ui - Vi 



\2H-2 



[0,1] J[u,u+1] 



M(dui)M(dvi) 
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(p-r)/2 



X 

i=l 



f \w-Vi\™- 2 M(dui)M(dvi) 



x n / \ui-vi? H - 2 

«=(p— r)/2+l 



Set 



x M(dui)M(dvi) 

e u (r,p)' := m 2pff-«iK2p) E / TT l , l , ^ 



ml 2 *" 2 



r+(p-r)/2 

3, is 

i=r+l 
P 

X l A a |tii-«i| 21r - 2 

i=r+(p— r)/2+l 

As in step 2, this expression is bounded and does not depend on n. Also, 
since there are r differences with — Vi\ > u/m n , there exists a constant C 
such that for all u |Jj"^"]fj_ 2) < C. 
Define 

e u (r,p)'m^ 1-a ^ 



I/n(«) 



r{2H-2) 



Under A(p), U n {u) converges a.s. to a certain random measure W U)T . := 
l^"(2H-2) x M^ 2p )([0, 1]). Moreover, as in step 2 of II, the random part of 
this measure does not depend on u. Also, define 

D n (u)=E Pn (Wo, n )-U n (u). 

In this case bounding the expectation of 

u-l 

E(\w ln (v)-w ln (j/m n )\) <Cj2HMU + k)/m n )-w ln U/m n )\) 

k=l 

+ E(\w ln (v)-w ln ((j + u-l)/m n )\) 
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will yield 

E(\D n (u)\) = 0( m a-°)(V-(2pr)-^(2 P ))/r-a(r-l)/r-^_ 

IV. To end the proof, we must show that m^ 3H ^^ 2p ^F n (r,p) converges 
a.s. to the limiting r.v. T(r,p) = W + 2Y^=iW Uir u r ^ H ~ 2 \ For this, it is 
enough to consider fas ir(2ff-2) is 

the term of a convergent seriesl 
Let A n be as in I. It is enough to show that 

P\\^ r(2H ' 2) \Wu,n ~ U n (u)\ > £ j n jCj 

is the term of a convergent series. This follows from remarking that, for any 
m, the above probability is bounded by 

C ' m (7-l+a)/2+(l-a)(^(2pr)-r^(2p))/r-a(r-l)/r m 2 

+ sup E(W u , r + ^,n)m r(2H ~ 2)+1 

and that E(W„, r ) < CE(W . r ) and E(W 1i) „) < CE(W 0>n ). 

The limit T(r,p) is integrable because Wq is and is a.s. positive as T(r,p) > 
W . □ 

Let g r be the coefficients of the expansion of G(x) = \x\ p — c p over the 
Hermite polynomials (since G is a centered even function, g r = for r = 0, 1). 
Set T(p) = J2 r gr^-T(r,p). Remark that since J2 r 9r r * < 00 anc ^ r n (r,p) < 
T n (2,p), we have r n (p) :=Erflr r '^n(''!P) converges a.s. to T(p), where the 
latter is a positive integrable r.v. 

Lemmas 3.1 and 3.2 assure the conditional convergence of ^( 2p )/ 2 Z n (p) 
toward a normal random variable. More precisely, we have the following: 

Theorem 3.1. If 1/2 < H < 3/4, then for even p, under Hl-2p and 
condition A(p), given in Lemma 3.2, then a.s. 

E M (e rf "* (2p)/2z "( ! '») e -^ 2r W/ 2 . 

Proof. Write 

pH-ifi(2p)/2 

m pH-^ P y2 Zn{p) = mn £ ^ G(AXfc/afc|n)) 



l<fc<m„ 



where, by Lemma 3.1, {AXk/ak,n}keN is a (conditionally) Gaussian and cen- 
tered sequence of strongly dependent r.v., that is, such that Cov(AXfc/afc ]n , 
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AXj/a jjn ) = 0(\k - j\ 2H ~ 2 ). The result follows directly since G has Her- 
mite rank equal to 2, [3, 9] by calculating the conditional moments of Z n (p) 
and by the a.s. convergence of T n (r,p). Indeed, using moment calculating 
techniques, it is possible to show that 

E M ([mV H -^/ 2 Z n (p)] 2 < 1 ) = c 2q T n {pf(l + o n (l)) 

and 

E M (Kf -^/ 2 Z n (p)] 2 « +1 ) < o n (l)T n {p)^ +l)/2 - 
The proof then follows by the a.s. convergence of T n (p). □ 

4. Statistical applications. In this section we study a moment based es- 
timator of C(p)- Moment estimators have the advantage of being very easy 
to implement. However, experimental studies show the existence of a bias 
problem as discussed in Section 1. 

Based on the results of Section 3, we will discuss rates of convergence for 
these estimators and offer insights as to some of the estimation problems at 
hand. Namely, we are interested in the construction of asymptotic confidence 
intervals. When m n 3> T (T = 1 in order to simplify), we show the existence 
of a random l/logm n bias term, whereas the conditional CLT yields an 
estimation error of order l/(ml/ 2+ ^^ ^ 2p ^ 2 log(m n )). 

When T ~ m n , the scaling property of the MRM yields the appearance 
of a linear bias term, as shall be discussed in detail in Section 4.2. In this 
case approximate confidence intervals are given. 

Based on the notation introduced in Section 3, T>m will stand for conver- 
gence in distribution conditional to the MRM M. 

4.1. Case T = 1. Throughout this section, we will assume p is an even 
positive integer. Set to begin with 

pH-4>{p) m n 

and 

P H -i/>(p) rn n 
Bn{p) = — C P a ln- 

We have a first result concerning the convergence of B n (p) . 



Lemma 4.1. Assume Hl-p holds. Then, B n (p) converges a.s. to a cer- 
tain a.s. positive random variable B(p). 
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Proof. Convergence of B n (p) to the random measure B(p) = e' (p/2) x 
M(p) [0, 1] follows as in the proof of Lemma 3.2. □ 

Using a first-order Taylor expansion and Theorem 3.1 yield, whenever 
if)(2p) < 2ip(p) + 1 

m l/2+^(p)-^(2p)/2 lQg f Snjp) \ = m l/a+V,(p)-V.(2p)/a S n -B n (p) 
\ B n J B(p) + o n (l) 

y M N(0,T(p)) 
B(p) 

The above discussion yields the following result. 

Corollary 4.1. Assume Hl-2j> holds true for p an even integer. Under 
the conditions of Theorem 3.1, 

log S n {p) , 



logm^ 

„^_l/2+^(p)-^(2p)/2 1 ^^ / _ , 



where log(m n )b n — > \og{B(p)) a.s. andml/ 2+ ^^ ^( 2p )/ 2 log(m TO ) „ B ^\ e r , - 
iV(0,l). 



4.2. Case T ~ m n : a certain linearization effect. In this subsection we 
assume that T = 0{m n ). That is, that the range of scale invariance of the 
MRM M is large enough to be comparable to the (fixed) number of obser- 
vations. In this case the scaling property accounts for a certain linearization 
phenomena in the estimation procedure. 

We are interested in nonasymptotic bounds for | log(5 n (p))/log(m n ) — 
E(log(5' ri (p))/log(m n ))|. More precisely, we shall see in this case that the 
scaling property yields E(| log(5 n (p))/log(m n ) — (H + i/j'(0))p\) < C j \f\ag~m r _ 
assuming ip is two times differentiable at zero. 

We must first introduce some additional notation. Define 

ri ri l 1 / 2 

\u-v\ 2H - 2 M{du)M(dv) 



j-Uj-l 
and 

SUlti\u-v\ 2H - 2 M{du)M{dv) 

PlJ, k) = — 

aja k 

By the scaling property, 4 in Lemma 2.1, we have 

p(j, k) l = p(k,j), 
ml H a kjn a jtn l = e 2Cl / m ™ d k dj. 
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Next introduce the random matrix S n := [a,kcijp(j,k)]j,k and consider S n 
the conditional covariance matrix of m„ (AXi, . . . , AX n ). The scaling prop- 
erty thus yields 

For each n, let (AXi, . . . , AX„) be a conditionally Gaussian centered ran- 
dom vector with conditional covariance matrix S n . By the scaling property 

m* (AX! ,...,AI n )^ w e Q v™« (AX a , . . . , AX n ) . 

Remark that, as follows from the scaling property, the same bounds hold 
a.s. for p as for p. Define for r S N, r > 2 



T n (r,p) = — ]T p{k,j) r a r -a p k . 



m ™ l<k,j<m n 

Then, by stationarity of the MRM, with the notation introduced in the 
proof of Lemma 3.2, 



E(m^f n (r,p)) < C r e (p) 1 + 2 £ 



u r(2H-2) 



U>1 



whenever Hl-2p holds. 
Introduce, for p£N, 

-i m n 

By the scaling property, 

Z n (p) ^ePni/mn^Cp), 

where, by construction, the random variable ^i/ mn is independent both 
of the collection of random variables {M[k,k + l|}feez for k < m n and the 
fractional Brownian motion . 

Calculating the conditional variance, 

E M ((Z n (p)) 2 )=]>> 2 r!f n (r,p), 

r 

so that under Hl-2p, we have 

E(m% H (Z n {p)f) < e (p) $> 2 r!^ ( 1 + 2 £ u^" 2 )) . 
Next consider, for even p, 
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For each fixed n, B n (p) is a bounded and positive r.v. under Hl-p, such 
that there exists a constant ki, E(B n (p)) < kieo(p/2). Jenssen's inequality 
and positivity of eg(p/2) thus yield that there exists a constant ki with 
|E(log(S n (p)))| <k 2 . Let 



pH 



5n(p) = ^El A ^i 



.IP 



™n ■ 'i 



Once again the scaling property assures 

I S n \ l aw / S'f) 

log( jd =log U 



On the other hand, 

log(J3„) = log(5 n ) +p^i /mn , 

where ^i/ m , n is an infinitely divisible random measure, independent of B n 
and of the fractional Brownian process Bf 1 , and such that ~E(e q ^ x ) = e — log ( A ) . 
Note that the latter yields 

with and independent r.v. Hence, if m = \ogm n , writing 
l/m n = ((l/m n ) 1/m ) H ((lMn) 1/m ) m - [ml , 

we have 

^l/mn L = ^'l H H ^[logm n ] + ^' 

where f2j are i.i.d. r.v. with E(e 9n j) = e^^\ Let q be such that > 1. By 
definition, E(e 9 "*) = e -V>(<7)(m-M)log(i/e)_ B y Markov's inequality, 

e ip{q)m 

P(n* > m) < = e-tt q)m , 

which is the term of a convergent series and f2*/logm n — > a.s. It follows 
([11], page 565) that 1 / mn /log(m n ) — ► ^'(0) in probability whenever ^'(O) 
exists. Moreover, if i/j"(0) exists, then E((f2') 2 ) < oo ([11], page 512), which 
by the SLLN assures U 1 / nin / log(m n ) — > tp'(0) a.s. and that 

n v™» ^/ (0) " 



m r , 



N(0,a 



2 



(9) v^og.. 

logm 

by the CLT, where a 2 n = Var(^i). 

An application of Markov's inequality and the above discussion yield the 
following result. 
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Corollary 4.2. Assume ip"(0) exists and set /j, := ip' (0) . Assume Hl- 
2p holds true. Then if m n = 0(T), there exists a positive constant C such 
that 
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